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Abstract. We present the results of an extensive timing analysis of the 1999 outburst of the soft X-ray transient 
and black hole candidate XTE J 1859+226 as observed with the Rossi X-Ray Timing Explorer. Three main different 
types of low frequency (1-9 Hz) quasi-periodic oscillations (QPOs) were observed and classified, strenghtening 
the general picture that is emerging for the variability properties of black hole X-ray binaires. Rapid transitions 
between different power spectral shapes were observed and their link with the count rate was studied. Furthermore, 
we show that a frequency of ~6 Hz seems to hold a particular place: one of the three QPO types we found was 
very stable when at this frequency, as it happens in Z sources as well. The coherence of its subharmonic peak was 
higher when the fundamental was close to 6 Hz, thus suggesting the presence of some resonance at this frequency. 
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' 1. Introduction Ivan der Klis 19951 for an overview). While only a few 

. . , detections were reported on the basis of EXOSAT data 

Observations with the Rossi X-ray limvnq Explorer . , on , , lr , , , , n>-nr\ 

, i , ,. • , m the early 80 s, the Gmga satellite showed clear QPO 

(RXlbj) have led to an extraordinary progress m the c , , ., mn n> c titt/-i 

, ; . . J rl , , , features m the power density spectra (PDb) of BHCs 

knowledge .of the variability properties o^ck-h^kcanc^ (gee | van der Klis 1995 1 . With RXTE, low-frequency 

dates (BHCs) m X-ray binaries (see e.g. van der Khs 2000 r\r>r\ j , j • n n u j n TT /~i 

lPr H — n — . . „ nnn i _., . — : : — ~ QPOs were detected m virtually all observed BHCs 

and Rcmillard et al. 2002a). 1 he fast quasi-periodic oscil- , i t — ttt- — norm tt n r\r,r\ 

. . 7T-— — , * , , (see [van der Khs 2004). Usually these QPOs are asso- 

lations QPOs that were discovered in many of these sys- . , , •,, 7 , , • , •>• , , , 

v . , r ciated with the spectrally hard and intermediate states 

terns are thought to originate in the innermost regions of i van der Klis 1995 1 iMcClint oerTO& lard 2004): they 

the accretion flows around stellar-mass black holes. Even ,, .,, a , . , 

. ^t-.^ ■ -ii appear together with a hat-top noise component, and 

though the mechanism responsible for the QPOs is still , . , r , , i. • <• a. 

, . . I,,. on a time scale of days their frequency often corre- 

unknown, the study of their properties and behaviour can , , ., , , , , / i /h ■ mno i j 

. . J ii. r • lates with the source count rate (see e.g. Cm f 998) and 

provide important clues on the physics of accretion onto rR— : t — -, „ ^ -u 

gjj^ , 1 * J |Keig et al. 2000| |. However, m the Gmga observations 

. of the bright transient GS 1124-683 two distinct types 

While only a few BHCs show high-frequency c , c r^nr\ -j i-c j • j. j 

J to H J of low frequency QPOs were identified: one associated 

QPOs (50-450 Hz, for a recent review see .,i n , , • , j , , 

■■T*. i m i 1. m - „ ' ■ . with a nat-top noise component and one to a steep noise 

McClmtock & Rcmillard 2004 and references therein), , , . j — i t f.^^ T , c , , , , 

■ — — 1 ; ' component ( lakizawa ct al. 1997). Ihe first type showed 

low- frequency QPOs with frequencies ranging from a few , , , , , , , , 

M J ^ H ° ° a strong dependence on count rate, while the other had 

mHz to ~10 Hz are a common feature. The low-frequency , £ T i , r^T~>r\ u j j. -j 

M J a rather stable frequency, lhese two QPOs had centroid 
QPOs were already known before the RXTE era (see r . 1 1f1 n , , n , , , 

_2 J v frequencies m the 1-10 Hz range and were clearly related 

Send offprint requests to: casella@mporzio.astro.it to the two PDS 'flavors' of very h igh state observed in 

* Present address: INAF - Osservatorio Astronomico di this system ( |JVliyamoto et al. 1993| ). 

Brera, via E. Bianchi 46, 1-23807 Merate (LC), Italy 
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Fig. 1. RXTE/PCA light curve (PCUs and 2, top panel), hardness ratio [(7-15 keV)/(2-7 keY),middle pane\ and 
total rms (0.03-64 Hz, bottom panel) of XTE J1859+226 during its 1999 outburst. The time resolution for the light- 
and hardness curves is 16 s except for the last five points, which have one point for each observation. 



Wijnands et al. (1999) and Homan et al. (2001) re- 
ported on two different types of QPOs in the RXTE data 
of XTE J1550-564: a broad one (type-A), with a quality 
factor Q (the QPO frequency divided by the QPO full- 
width-half- maximum (FWHM)) of less than 3, and a nar- 
rower one (type-B), with a Q larger than 6. Both QPOs 
were characterized by a centroid frequency of 6 Hz and as- 
sociated with a weak red- noise component, but with differ- 
ent phase-lag behaviours. XTE J1550-564 also showed the 
more common QPO-type associated with a flat-top noise 
component (see ICui et al. 19991 and ISobczak et al. 20 00 ) . 
IRemillard et al. 2002bl dubbed this QPO 'type-C: its fea- 
tures are a high coherence (Q 10), a variable centroid 
frequency (in the range 0.1 - 10 Hz) and a strong flat-top 
noise component (~10-40% rms) (see Table 0. 

While type-C QPOs are observed in many systems, the 
other two are less common. In addition to XTE J 1550- 
564, type-B QPOs were also observed in GX 339- 
4 (e.g. |Miyamoto et al. 1991| and |Nespoli et al. 2003| |, 
GRS 1739-278 (|Wijnands et al. 2001|) and possibly in 
4U 1630-47 (|Tomsick fc Kaaret 2000"|l, while type-A 
QPOs were observed in GX 339-4 flNespoli et al. 2003|> 



and possibly in 4U 1630-47 l|Tomsick & Kaaret 2000f aud 
IDieters et al. 2000)l . Furthermore, in the light of the A-B- 
C classification, the two QPOs observed in GS 1124-683 



Table 1. S 



ummary of type-A, -B and -C QPOs p rop erties in XTE J1550-5 64 
)Wijnands et al. 1999llHoman et al. 20011 and IRemillard et al. 2002bl 



Property 



TypeA 



TypeB 



TypeC 



Frequency {Hz) 
Q {v I FWHM) 
Amplitude (%rms) 
Noise 



~ 6 
<3 
3-4 
weak red 



~ 6 
>6 

4 

weak red 



0.1 - 10 
> 10 
3-16 
strong flat-top 



(see above) can be tentatively identified with types B and 
C, although a detailed analysis of Ginga data is necessary 
to confirm this association. 

These oscillations, whose nature is still not under- 
stood, provide a direct way to explore the accretion flow 
around black holes (and neutron stars). In particular, 
their association with specific spectral states and the phe- 
nomenology that is emerging indicate that they are a key 
ingredient in understanding the physical conditions that 
give origin to the different states. 

The soft X-ray transient XTE J1859+226 was discov- 
ered on 1999 Octo ber 9 with the RXTE All Sky Monitor 
<|Wood et al. 1999jl . which detected a 2-12 X-ray flux of 
^160 mCrab (~ 4.5 x 10 -9 erg s _1 cm' 2 ) quickly ris- 
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Table 2. Power-spectral classification and variability parameters' 1 . 
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Obs. Id. 


Date 


MJD 


QPO Type 




FWHM qpo 




PhaseLag c 


0.03-64 Hz rms b 












(Hz) 


(Hz) 


% 


(rad) 


(%) 


40124-01-04-00 


1999 Oct 11 


51462.768 




c 


1.2 


0.18 


16.61 


0.00±0.01 


29.8±0.3 


40124-01-05-00 


1999 Oct 12 


51463.833 




c 


3.05 


0.31 


13.48 


-0.04±0.01 


24.9±0.2 


40124-01-06-00 


1999 Oct 13 


51464.109 




c 


3.64 


0.46 


13.43 


-0.07±0.01 


22.9±0.1 


40124-01-07-00 




51464.633 




c 


3.65 


0.46 


12.97 


-0.04±0.01 


23.1±0.2 


40124-01-08-00 


1999 Oct 14 


51465.305 




c 


4.39 


0.48 


11.85 


-0.08±0.01 


20.78±0.08 


40124-01-09-00 




51465.498 




c 


4.95 


0.56 


10.85 


-0.09±0.01 


19.16±0.08 


40124-01-10-00 




51465.902 




c 


5.81 


0.72 


8 


.52 


-0.14±0.01 


16.2±0.1 


40124-01-11-00 


1999 Oct 15 


51466.896 




c 


5.97 


0.65 


8 


11 


-0.18±0.01 


16.0±0.08 


40124-01-12-00 


1999 Oct 16 


51467.581 




A 


7.75 


3.85 


1 


.75 


-0.33±0.07 


2.76±0.05 


40124-01-13-00 




51467.961 




B 


6.08 


0.45 


2 


.10 


0.07±0.07 


4.5±0.06 


40124-01-13-00 




51467.973 




A 


8.42 


2.885 


1, 


.62 


-0.3±0.1 


2.8±0.1 


40124-01-14-00 


1999 Oct 17 


51468.427 




C* 


8.64 


2.37 


2 


.83 


-0.56±0.04 


7.67±0.04 


40124-01-14-00 




51468.427 




B 


6.43 


0.82 


3 


13 


0.3±0.09 


6.1±0.1 


40122-01-01-03 


1999 Oct 18 


51469.093 




B 


5.96 


0.75 


3 


04 


0.23±0.04 


5.02±0.04 


40122-01-01-02 




51469.159 




B 


5.96 


0.85 


2 


52 


0.13±0.04 


4.83±0.05 


40122-01-01-02 




51469.179 




A 


7.85 


3.26 


1 


.70 


-0.02±0.1 


2.9±0.1 


40122-01-01-01 




51469.226 




A 


7.74 


5.22 


1 


.37 


-0.84±0.2 


2.1±0.1 


40122-01-01-00 




51469.293 




A 


7.59 


3.81 


1 


64 


-0.33±0.1 


2.16±0.09 


40122-01-01-00 




51469.360 




B 


5.94 


0.80 


2 


.02 


0.125±0.06 


4.40±0.06 


40124-01-16-00 




51469.493 




C 


8.45 


3.54 


3 


10 


-0.49±0.04 


7.69±0.07 


40124-01-16-00 




51469.517 




c* 


7.7 


1.36 


4 


.53 


-0.29±0.04 


8.28±0.09 


40124-01-17-00 




51469.893 




c* 


7.89 


1.17 


3 


.22 


-0.41±0.07 


9.6±0.1 


40124-01-17-00 




51469.899 




c* 


7.54 


1.58 


5. 


07 


-0.26±0.04 


10.6±0.1 


40124-01-17-00 




51469.905 




c* 


7.27 


1.19 


5. 


.40 


-0.22±0.04 


11.7±0.1 


40124-01-18-00 


1999 Oct 20 


51471.024 




c 


5.88 


0.42 


6 


98 


-0.17±0.015 


15.81±0.07 


40124-01-19-00 




51471.224 




c 


5.18 


0.61 


9 


10 


-0.15±0.015 


18.11±0.07 


40124-01-20-00 




51471.890 




c 


6.47 


0.80 


6 


.59 


-0.25±0.025 


14.04±0.07 


40124-01-21-00 


1999 Oct 21 


51472.503 




c 


6.34 


0.90 


7 


20 


-0.2±0.02 


14.49±0.06 


40124-01-15-00 


1999 Oct 22 


51473.245 




c 


6.13 


0.72 


7, 


11 


-0.22±0.04 


15.5±0.1 


40124-01-23-00 




51473.822 




c 


6.92 


0.93 


5. 


.68 


-0.27±0.03 


12.79±0.09 


40124-01-23-01 




51473.890 




c* 


7.73 


1.59 


4 


.00 


-0.43±0.03 


9.85±0.07 


40124-01-15-02 


1999 Oct 23 


51474.088 




c* 


7.39 


4.52 


3 


63 


-0.43±0.03 


7.91±0.08 


40124-01-24-00 




51474.429 


B 


— Cath. 


5.84 


0.81 


3 


.62 


0.23±0.04 


6.27±0.06 


40124-01-25-00 




51474.820 




c 


6.19 


0.68 


7, 


.18 


-0.22±0.02 


14.9±0.1 


40124-01-26-00 


1999 Oct 24 


51475.154 




C* 


7.06 


0.95 


5 


.06 


-0.29±0.03 


12.28±0.08 


40124-01-26-00 




51475.166 




c* 


7.24 


1.06 


4 


.43 


-0.34±0.03 


11.31±0.07 


40124-01-26-00 




51475.177 




c* 


7.73 


1.36 


3 


.54 


-0.40±0.03 


9.84±0.07 


40124-01-26-00 




51475.218 




c* 


7.75 


1.72 


3 


.70 


-0.45±0.02 


9.29±0.04 


40124-01-26-00 




51475.252 




c* 


7.56 


1.24 


3 


.65 


-0.33±0.03 


10.35±0.07 


40124-01-27-00 




51475.428 


B 


— Cath. 


5.79 


0.84 


4 


.28 


0.25±0.07 


6.2±0.1 


40124-01-28-00 


1999 Oct 25 


51476.428 




C* 


7.64 


2.1 


2 


.91 


-0.48±0.03 


8.69±0.06 


40124-01-28-01 




51476.501 




c* 


7.52 


1.53 


4 


.22 


-0.36±0.04 


10.27±0.08 


40124-01-29-00 


1999 Oct 26 


51477.152 




c* 


7.54 


2.4 


3 


.66 


-0.42±0.02 


8.99±0.04 


40124-01-30-00 


1999 Oct 27 


51478.017 




B 


5.06 


0.7 


3 


.70 


0.17±0.04 


5.36±0.05 


40124-01-31-00 




51478.777 




c* 


7.26 


0.96 


3 


.97 


-0.33±0.03 


11.12±0.07 


40122-01-02-00 




51478.816 




c* 


7.77 


1.31 


2 


74 


-0.40±0.05 


9.51±0.08 


40122-01-02-00 




51478.844 




c* 


7.77 


1.89 


3 


.38 


-0.44±0.05 


8.9±0.1 


40122-01-02-00 




51478.882 




c* 


7.38 


1.25 


4 


.00 


-0.36±0.04 


10.72±0.09 


40122-01-02-00 




51478.913 




c* 


7.54 


1.39 


3 


.59 


-0.40±0.04 


10.45±0.08 


40122-01-02-00 




51478.949 




c* 


7.04 


0.96 


4 


84 


— 0.37±0.03 


11.99±0.07 


40122-01-02-00 




51478.980 




c* 


7.46 


1.03 


3 


12 


-0.49±0.05 


10.01±0.08 


40124-01-36-00 


1999 Nov 01 


51483.117 




B 


4.69 


0.63 


4 


.00 


0.16±0.03 


5.58±0.07 


40124-01-36-00 




51483.158 




D 


4.69 


0.54 


3 


93 


0.22±0.04 


5.88±0.06 


40124-01-37-00 




51483.945 




B 


4.59 


0.54 


4 


31 


0.16±0.03 


5.80±0.05 


40124-01-37-01 


1999 Nov 02 


51484.077 




B 


4.53 


0.62 


4 


.05 


0.07±0.05 


4.89±0.15 


40124-01-37-02 




51484.276 




B 


4.45 


0.61 


4 


.03 


0.09±0.03 


4.33±0.06 


40124-01-39-00 




51485.875 




? 


5.61 


2.3 


2 


.68 


-0.51±0.04 


4.1±0.1 


40124-01-40-00 


1999 Nov 04 


51486.828 




? 


5.37 


1.7 


3 


.00 


-0.43±0.07 


3.1±0.3 


40124-01-40-01 




51486.873 




? 


5.22 


1.3 


2 


.80 


-0.3±0.1 


4.0±0.2 


40124-01-41-00 


1999 Nov 05 


51487.009 




? 


4.74 


1.3 


2 


18 


-0.42±0.05 


3.5±0.1 


40124-01-42-00 


1999 Nov 06 


51488.409 




? 


7.08 


2.4 


1. 


.75 


-0.88±0.06 


4.46±0.08 


40124-01-49-01 


1999 Nov 15 


51497.255 




? 


5.21 


1.9 


2 


.82 


-0.36±0.07 


4.1±0.2 


40124-01-58-01 


1999 Dec 15 


51527.004 




? 


9.1 


0.8 


2 


17 


0.0±0.2 


8.9±0.2 


40124-01-59-00 


1999 Dec 18 


51530.133 




? 


8.15 


2.1 


2 


.91 


0.0±0.2 


6.71±0.25 


40124-01-61-00 


2000 Jan 03 


51546.036 




? 


4.65 


0.8 


2.3 


-0.17±0.25 


4.3±0.4 



a Only observations with evi dence for low frequency Q PO are listed. 
b Normalization according to Bclloni & Hasingcr 1990 

c Phase lag in radians between the 2-5 and 5-13 keV light curves and integrated over the frequency range u qvo ± FWHM/2 



ing at a rate of ~6 mCrab/hour. A follow-up observa- 
tion of the source with the RXTE/PCA (Proportional 
Counter Array, 2-60 keV) revealed a hard power-law dom- 
inated spectrum l|Markwardt et al. 1999|l . On October 16, 



the source reached its peak flux of ~ 5 x 10 -8 erg s _1 
cm~ 2 in the 2-80 keV band (corresponding to a lumi- 
nosity of ~ 7 x 10 38 erg s _1 for an assumed distance 
of 11 kpc (Zurita et al. 2002)). After the initial hard 
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Frequency (Hz) 



Fig. 2. Examples of type A, B and C QPOs from our 
XTE J1859+226 observations. QPO and harmonics peaks 
are indicated. Upper panel: obs. 40124-01-12-00. Middle 
panel: obs. 40122-01-01-03. Bottom panel: obs. 40124-01- 
10-00). The Poisson noise was not subtracted. 



rise, XTE J1859+226 softened at its peak intensity and 
continued to soften for almost two months, when a sec- 
ondary hard plateau took place (see Fig.^ MJD^ 51520- 
51540). This was rather similar to the behaviour observed 
in many X-ray transients, in particular XTE J1550-564 
(see IRemillard et al. 2002bl and references therein), sug- 
gesting a common scenario for the evolution of these ob- 
jects. Both low-frequency (~ 1 — 4 Hz and ~ 6 Hz) and 
high-frequency (~ 150 — 187 Hz) QPOs have been re- 
ported (see ICui et al. 20001 and IFocke et al."2000jl from 
XTE J1859+226. 

Th e optical jH arnavicn, Stanek, fc Berlind 1999) ) and 
radio ( |Pooley ^B^e^mm^^9^ counterparts were iden- 
tificd soon after the discovery of the source. Optical 
photometry revealed a possible period of 9.15 hr 
( |Garnavich &: Quinn 2 000). Monitoring at radio wave- 
lengths suggested that at least two relativistic ejec- 
tion episodes took place approximately on October 16.5 
(MJD = 51467.5) and 27 (MJD = 51478). However, 
no ejecta were spatially resolved at radio wavelengths 
HBrocksopp et al. 2002| ). 

Here we present the results of an extensive X-ray tim- 
ing analysis of the 1999 outburst of XTE J1859+226, fo- 
cussing on the low-frequency QPOs. We found three dif- 
ferent types of low frequency (1-10 Hz) QPOs. We show 
that these correspond to the above-mentioned A, B and 
C QPOs types, all of them showing distinctive and well 
defined behaviours and phase lags . 



Table 3. RXTE/PCA data modes active during the 
XTE J1859+226 observations 



Mode 


Time res. 


Number of 


PHA 


Name 


(sec.) 


PHA Channels 


Energy range (keV) 


Standardl 


2- 3 


1 


2-60 


Standard2 


2 4 


129 


2-60 


Binned 


2- 7 


36 


2-15 


SB1 


2 -13 


14 


2-6 


SB2 


2 -13 


22 


6-15 


Event 


2 -16 


16 


15 - 60 



2. Observations and data analysis 

We analyzed 129 RXTE/PCA observations made dur- 
ing the 1999 outburst of the black-hole candidate 
XTE J1859+226, between MJD 51462 (1999-10-11) and 
51626 (2000-03-23). Table |2 shows dates and parameters 
of the observations where a low frequency QPO has been 
observed. 

The PCA data were obtained in several simulta- 
neous different modes (see Tab. El • Only proportional 
counter units (PCUs) and 2 were always active dur- 
ing our observations. Standard 2 data from these two 
PCUs were used to create light- and hardness curves 
for the whole outburst, whereas the high time resolution 
data from all active PCUs (in a given obseravtion) were 
used for the timing analysis. A hardness ratio was de- 
fined as the ratio of counts in the range 7-15 keV (12- 
31 channels) to those in the range 2-7 keV (0-11). Fast 
Fourier Transforms were made from 16s data intervals 
with Nyquiest frequencies of 64 Hz (Binned data) and 
4096 Hz (Single bit data). The resulting PDS were 
averaged, rebinned logarithmically, and the Poissonian 
noise, including the Very Large Events (VLE) contri- 
bution jZhang 1995| |Zhang et al. 1995| , was subtracted. 
The PDS were normalised to fractional squared rms, 
following |Belloni fc Hasinger 1990| PDS fitting was car- 
ried out by using the standard Xspec fitting package 
by using a one-to-one energy-frequency conversion and a 
unit response. Following Bcllo nTet al. 20021 we fitted the 
noise components with two Lorentzian shapes, one zero- 
centered and a second one centered at a few Hz. The QPOs 
were fitted with one Lorentzian each too, only occasionally 
needing the addition of a Gaussian component to better 
approximate the shape of the narrow peaks and to reach 
values of reduced x 2 close to 1. For the observations where 
the dynamical power spectra showed transitions between 
different power spectral shapes (see below), we separated 
different time intervals in order to obtain average power 
spectra for each shape. 

For every 16 s interval we also produced a cross- 
spectrum between the 2-5 and 5-13 keV resolved light 
curves (defined as C(j) = X^(j)X 2 {j), where X\ and 
X2 are the complex Fourier coefficients for the two en- 
ergy bands at a frequency i/j), calculated average cross- 
spectrum vectors for each observation, and then derived a 
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Fig. 4. Same as Figure Q but for only the first 25 days 
of the outburst. See inset for the symbols that are used 
to represent different QPO types in the lower panel - dots 
are used if no QPO was detected. 



phase lag as a function of frequency from the angle in the 
complex plane of these vectors {<pj = arg[C(J)]). The error 
in <f> is computed from the observed variance of C in the 
real and imaginary directions. In line with recent litera- 
ture we defined phase lags as positive when the hard X-ray 
variability follows the soft one. To quantify the phase-lag 
behaviour of the QPOs, we extracted the phase lags in 
a range centered at the QPO peak frequency and corre- 
sponding to the width of the peak itself (y v ± FWHM/2, 
see |Reig et al. 2001)1 ) . In Tabled we list the frequency, full- 
width half- maximum (FWHM), 2-15 keV fractional rms, 
and phase lag of the QPO for each observation in which 
one was detected. The total integrated fractional rms (2- 
15 keV) of the PDS is given as well. 

In Figure m we show the 2-60 keV light curve of the 
whole outburst, the hardness ratio and the integrated 
0.03-64 Hz fractional rms of the 2-15 keV light curves. 

3. Power spectra and phase lags 

In many of our observations QPOs were detected, with 
frequencies ranging from ~1 to ~9 Hz. Three main types 
could be distinguished, which, based on their phase lag 
and coherence properties, could be associated with type- 
A, -B, and -C QPOs. Example power spectra of each type 
are shown in Figure |21 In addition to these three main 
types we also identified various sub-types, which will be 
discussed below. 

A useful method for differentiating between the three 
types of PDS is shown in Figure [3J in which we plot 



the integrated fractional rms of each PDS versus the cen- 
troid frequency of the QPO. Several groups of points can 
be identified. The first large group of points in the plot, 
type-C and C* (see paragraphs $3.1 and $3.2), is diago- 
nally spread across the plot, covering the whole frequency 
range between 1 and 10 Hz and a large range in rms (7- 
30% rms). Another group, type-A (see paragraphs $3.3), 
is clustered around a frequency of 8 Hz and rms of 2%. 
Finally a third group is located at a slightly higher rms (4- 
6%) in the 4.5 - 6.5 Hz range: type-B and " B-Cathedral" 
(see paragraphs $3.4 and $3.5). A more detailed analysis 
of the correlation between the frequency these three QPO 
types and the integrated fractional rms in different sources 
is described in a forthcoming paper (Casella et al. 2004, 
in prep.). In Figure 0] we show the light curve, the hard- 
ness and the total fractional rms of the first 25 days of 
the outburst (see Fig ^ for energy and frequency ranges) 
and indicate where the three types of QPO appear. In the 
bottom panel we have marked the three types of QPO and 
their associated sub-types. In the following paragraphs we 
describe all these types in detail, characterizing their PDS 
and phase-lag behaviour. 

3.1. Type-C QPOs 

In the early stages of the outburst (see Table [5] for dates), 
when the source flux was rising quickly, the PDS showed 
four main components: a strong (15-30% rms) flat-topped 
noise and three harmonically related QPOs (see Fig[3J). 
The central 1-7 Hz QPO was strong (rms amplitude 6- 
16%) and narrow (Q ~7-10). Note that in a few cases 
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Fig. 5. An example of type-C power spectrum (2-15 keV; 
obs.: 40124-01-10-00). The solid line shows the best fit 
with five Lorentzians (dotted lines). See the correspondent 
lags in Fig. [U fourth panel. 



the addition of a Gaussian component was required, in 
order to better approximate the peak's shape. This was 
also the case for the second harmonic, which was always 
present. A broad (Q ~2) subharmonic was also detected 
in all cases, except for the first observation (40124-01-04- 
00) where the frequency of the fundamental was at its 
lowest value). The harmonic relation between the three 
peaks was confirmed by allowing the centroids to vary 
independently. We decided to fix the centroid frequencies 
to harmonic ratios to obtain self-consistent estimates of 
the amplitude and width. 

The phase-lag behaviour of the type-C power spectra 
is similar to that found in GRS 1915+105 by Reig et al. 
(2000). It is strongly correlated with the frequency of the 
QPO, with a trend towards negative lags for increasing 
QPO frequency. Figure HO (all panels except the bottom 
one) shows four examples covering the whole observed 
range in QPO frequency. Owing to poor statistics, lags 
became unmeasurable at high frequencies; we thus plot 
them only below 20 Hz. In all observed cases, the funda- 
mental of the QPO showed negative lags (see Figure IT71 
and the Discussion), with a clear trend towards zero for 
decreasing centroid frequency, consistent with the Reig et 
al. (2000) results on GRS 1915+105 (where peak lags are 
positive at frequencies below ~1 Hz). The lags of subhar- 
monic peak were always negative as well, while the second 
harmonic always showed positive lags. It is important to 
mention the absence in the cross spectra of narrow features 
standing out at the frequency of the QPO and harmonics 
peaks. This suggests that we are measuring the phase-lags 
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Fig. 6. Type-C phase lags vs. frequency for four different 
QPO frequencies. Bottom panel: an example of type-C* 
phase lags. Positive values indicate that the hard (5-13 
keV) photons are lagging the soft(2-5keV) photons. QPO 
centroid frequency and observation I.D. are indicated for 
each panel. The dashed lines mark the frequency of the 
QPO, while the dotted lines mark the subharmonic (if 
present) and second harmonic frequencies. 



of the underlying noise continuum rather than the phase 
lags associated with the QPO peaks. 
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Fig. 7. An example of type C* PDS (2-15 keV) and phase 
lags (obs.: 40124-01-28-00). The dotted vertical line marks 
the frequency of the QPO. 




Frequency (Hz) 



Fig. 8. An example of type-A PDS (2-15 keV) and phase 
lags (obs.: 40124-01-12-00). The dotted vertical line marks 
the frequency of the QPO. Lags at frequencies above ~ 
10 Hz are, given the lack of statistics, random scattering 
between — ir and w. 

3.2. Type-C QPOs 

Later in the outburst, type-C QPOs appeared again (see 
TablcEJ. In this case, however, the centroid frequency was 
higher and the rms amplitude lower. We refer to these 
QPOs as a type-C*. The power spectra showed a strong 



red-noise component and a broad QPO peak at a cen- 
troid frequency of 7-9 Hz (Figured upper panel). Again 
a Gaussian component was added in many cases in order 
to better approximate the peak shape. A second harmonic 
peak was sometimes present, as well a subharmonic one. 
Phase lags were negative and large, up to 10 ms, over the 
range v p ± FWHM/2, and decreased rapidly at frequencies 
slightly higher than the QPO centroid, as can be seen in 
Figure (bottom panel). 

Even though the QPO centroid frequency range was 
different from that of type-C QPOs, and the rms and Q- 
values were smaller, the results above (see also Figure 
provide evidence that the properties of type-C and type- 
C* QPOs are smoothly connected if ordered for increasing 
QPO frequency. The phase-lag behaviour confirms this: 
even though the shape of a "normal" type-C QPO (Figure 
El upper panel) was rather different from that of type C* 
QPOs (FigurcHiJ bottom panel), a continuos transition be- 
tween the two took place for increasing QPO frequency. 
However, owing to the long time scale of the QPO fre- 
quency variability (not detectable during one single RXTE 
observation) it was not possible to observe a direct tran- 
sition between the two QPO types. 

3.3. Type-A QPOs 

Type-A power spectra (see an example in Fig. |S1 upper 
panel) appeared at the peak of the outburst (see Table 
|2J, when the count rate was very high. They were charac- 
terized by a broad QPO (Q ~ 2), with centroid frequen- 
cies between 7.5 and 8.5 Hz and fractional rms around 
~1.5%, and a low amplitude (few % rms) red- noise com- 
ponent. Neither a subharmonic nor a second harmonic was 
present. This was the QPO type with the lowest total rms 
in our sample. The phase lags, except for a negative excess 
around the QPO frequency, were consistent with zero. 

3.4. Type-B QPOs 

The narrow (Q ~10) QPO characterizing this type of PDS 
appeared only in a rather restricted frequency range be- 
tween 4.5 and 6.5 Hz (see an example in Fig. [5] upper 
panel). The red-noise was very weak (few % rms). The 
QPO peak profile was often more similar to a Gaussian 
than a Lorentzian, but we needed to combine both compo- 
nents in order to obtain a good fit (it is worth noticing that 
a similar combination was used by Homan et al. (2001) 
and Wijnands et al. (1999) for the XTE J1550-564 data, 
and by Nespoli et al. (2003) for the GX 339-4 data. In the 
latter case the authors explained the Gaussian shape with 
the presence of centroid frequency variability on a ~10 
s time scale. However, such variability was not detected 
in the XTE J1859+226 data). A weak second harmonic 
was always present, while a subharmonic appeared (with 
rms amplitudes <1 %) only when the fundamental fre- 
quency was at its highest values and the fundamental and 
the second harmonic had low amplitudes. This behaviour 
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Fig. 9. An example of type-B PDS (2-15 keV) and phase 
lags (obs.: 40122-01-01-03). The dashed vertical line marks 
the frequency of the QPO, while the dotted lines mark the 
subharmonic and second harmonic frequencies. 

suggests the presence of a sort of "balance" between the 
amplitude of the three peaks, particularly between the sec- 
ond harmonic and subharmonic: when the subharmonic 
peak was at its highest rms the second harmonic is not 
present, and viceversa. Unfortunately, poor statistics did 
not allow a more precise assessment of this. The lags (see 
Fig- IHl bottom panel) were positive at the frequency of the 
fundamental peak (see Fig. ^] for values) and negative at 
the frequencies of subharmonic and second harmonic. 

3.5. Type B- 'Cathedral' QPOs 

During two observations, 40124-01-24-00 (MJD: 
51474.429) and 40124-01-27-00 (MJD: 51475.428), the 
power spectrum showed a peculiar double-peaked QPO 
(see Fig ^3 upper panel) , similar to some of the type-B 
QPOs observed in X TE J1550-564 ( |Wijnands et al. 1999| 
Ho man et al. 200T|l . In both observations, two strong 
and narrow peaks were observed at harmonically related 
frequencies of ^3 and 6 Hz. In the case of the second 
observation we needed to add a Gaussian component 
to both peaks in order to better approximate their 
shape. A weak red-noise component was present (the two 
Lorentzians having ^2% rms each). The rms amplitudes 
of the two QPO peaks were ~4% and 2-2.5%, for the 
6 Hz and 3 Hz QPO, respectively. A weak (~1% rms) 
peak at the harmonically related frequency of 12 Hz 
was observed in both observations. The phase lags were 
consistent with zero over the whole frequency range 
except at the frequencies where the QPOs were seen (3 
and 6 Hz). The lags of the 3 Hz QPO were negative and 




0.1 1 10 

Frequency (Hz) 



Fig. 10. An example of type B-'Cathedral' PDS (2-15 
keV) and phase lags (obs.: 40124-01-24-00). The dotted 
vertical line marks the frequency of the QPO. 




Fig. 11. Coherence of the subharmonic peak for different 
pds types. 

corresponded (calculated, according with our definition, 
in the range v p ± FWHM/2) to a delay of ~20 ms in the 
soft X-ray variations, while for the ~6 Hz peak the lags 
were positive and corresponded to a delay of ^7 ms in 
the hard X-ray variations. 

One might classify the ~3 Hz peak as the fundamental 
and the ~6 Hz peak as the second harmonic. In this case, 
the ~12 Hz would be the fourth harmonic, while the third 
harmonic would be missing or too weak to be detected. 
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However, a comparison of this "cathedral" power spec- 
trum with the others, suggests a similarity with the type- 
B power spectrum. Apart from the high amplitude and 
narrower width of the ^3 Hz peak, the general character- 
istics are almost identical: the total rms was around a few 
% (compared to ~15% for the type-C spectrum when the 
QPO was at ~6 Hz), and both the red- noise and the ~6 
Hz peak were weak. Moreover, type-B PDS showed posi- 
tive lags for the fundamental, as opposed to those in the 
type-C power spectra. Thus it seems natural to identify 
the ~6 Hz peak (having positive lags) as the fundamental, 
the ~12 Hz as its second harmonic, and the ~3 Hz peak as 
its subharmonic. This is also in line with other PDS types 
in which we identified subharmonic peaks. Moreover, the 
integrated rms of the subharmonic peak of the B-cathedral 
PDS (~ 2 — 3%) was not much higher than that of other 
PDS. On the contrary, it was lower than that of type- 
C QPOs (between 3 and 8 %), and only slightly higher 
than that of "normal" type-B QPOs (~ 1 - 2%). However 
the coherence of the B-cathedral subharmonic peak was 
the largest observed among all low frequency QPOs (see 
Figure lTT|l . making the feature more prominent. It is worth 
noticing that the ~6 Hz QPO frequency clearly stands out 
in figure ITT1 fsee the discussion). 



3.6. Special cases: unclassified QPOs 

In a few of the observations following the last type-B PDS 
(MJD 51484) we detected QPOs that we were not able to 
classify in terms of the type A/B/C scheme. Owing to poor 
statistics and a too high centroid frequency variability, 
we could not unambiguously fit the power spectra and 
then obtain characteristic parameters. However, for sake 
of completeness we report in Table El the best estimate for 
the parameters of each observation. In Figure IT2*1 we show 
three representative examples of power spectra and phase 
lags from these observations. 



4. Timing and spectral evolution 

In some observations, the dynamical PDS showed rapid 
(a few tens of seconds) transitions between different power 
spectral shapes. In all cases, the transitions involve type-B 
QPOs. In panels a-c of Figure IT3l we show three examples 
of different behaviours. In the first half of the observation 
40124-01-13-00 (MJD 51467.961, panel a), when the light 
curve was highly variable, the PDS was of type-B (with 
a QPO frequency ~ 6 Hz). Simultaneously with the rise 
observed in the light curve after ~1100 s from the start, 
the PDS showed a sharp transition to a type-A shape with 
a QPO frequency of ~ 8 Hz (not visible in the gray scale 
representation). In the second part of the observation, the 
light curve was much less variable and had a higher mean 
count rate. Notice that a brief interval with the same char- 
acteristics (type-A QPO, higher flux) was seen ^200 s into 
the light curve, again with very sharp in and out transi- 
tions. 

For observation 40124-01-14-00 (MJD 51468.427, panel 
b) the behaviour was different: when the source flux was 
low the power spectrum showed a type-C shape (~ 8.7 
Hz), while during the two peaks in count-rate, when it 
reached values close to those of the first half of the previ- 
ous observation, the PDS was of type-B (~ 6.4 Hz). The 
transitions were again very sharp. In observation 40122- 
01-01-00 (MJD 51469.360, panel c) the source showed the 
opposite behaviour: at high count rates, the power spec- 
trum transitioned to type-A (~ 7.6 Hz), while at lower 
fluxes a type-B PDS was observed (~ 6 Hz), similar to 
panel (a). This was clearly seen before and after the gap 
in the middle of the observation. 

It is worth remarking that in all cases the transitions 
involve type-B QPOs. When the source showed a type-B 
PDS and underwent a fast transition to a lower count 
rate, the PDS changed to type-C; when on the other 
hand the transition was to a higher count rate, the PDS 
changed to type-A. In the same way, fast transitions from 
type-A to type-B and from type-C to type-B always in- 
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Fig. 13. Light curves and dynamical PDS for observations 40124-01-13-00 (panel (a), all five PCUs on), 40124-01-14-00 
(panel (b), five PCUs on), 40122-01-01-00 (panel (c), three PCUs on) and 40124-01-24-00 (panel (d), three PCUs on). 
The lower power value of the type-A QPO peaks with respect to the type-B render them invisible in the dynamical 
power spectra of panels (a) and (c). However, these are clearly seen in the total average PDS. 



volved a decrease and increase in count rate, respectively. 
Direct transitions between types A and C were not ob- 
served. Panel (d) of Figure [T31 (Obs. Id. 40124-01-24-00, 
MJD 51474.429), shows a fourth type of rapid transition, 
which occurred when the source showed a B-Cathedral 
type PDS. Corresponding to the dips in the light curve, 
the PDS changed its shape, with the two peaks par- 
tially losing their coherence while the red noise increased. 
Unfortunately, the time intervals in which this happened 
were too short for a detailed power-spectrum analysis. 



In Figure ITll we plot a light curve of all observations 
in which fast transitions were observed. From this figure 
it is evident that the count rate at which the transitions 
occur follow an exponential trend. Type-B QPOs appear 
in a narrow count rate range, as can be seen in the inset 
of the same figure, where we show the entire portion of 
the outburst where low frequency QPOs were detected. 

Similar fast transitions between different types of 
QPO and broad-band noise components were reported 
with Ginga from GS 1124-68 IjTakizawa et al. 1997jl and 
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Fig. 14. Light curve of the observations showing transi- 
tions between two different QPO types. The two paral- 
lel lines have been drawn so as to intersect the transi- 
tion points. In the inset the lines have been extended to 
the part of the outburst where QPOs appear. Different 
grayscales indicate different QPO types. (A color image 
is available on line) Black points correspond to data in 
which no QPOs were found. 



GX 339-4 itself flMiyamoto et al. 199l| ). Interestingly, also 
in these cases the sharp QPO has a frequency around ~6 
Hz. 

Finally, in Figure H3 (Obs. Id. 40124-01-38-01), we 
show a peculiar event: the source count rate dropped by 
^20% (with an integrated fractional rms of ~2.5% before 
and ~4.0% after the drop) in about 100 s. The increase 
in rms occurs mainly above 5 Hz, with a broad power ex- 
cess clearly visible around ~6 - 7 Hz (not shown). The 
hardness distribution was also different, as can be seen in 
the bottom panel of Figure El During the observation, 
the hard and soft colors decreased fairly continuosly, but 
simultaneously with the drop in count rate the soft color 
showed a clear discontinuity. This sharp transition is also 
clearly visible in Figuresn(MJD: 51485.075) where it can 
be seen that after the drop the count rate kept following 
an approximate exponential trend for several days, with 
a slope close to that shown in Figure 0] before the drop, 
which was thus a conspicuous feature in the overall light 
curve of the outburst. 

5. Discussion 

We analyzed the RXTE/PCA data from the 1999 out- 
burst of the black hole candidate XTE J1859+226, 
studying the low frequency QPOs and their detailed 
behaviour. We could classify most of QPOs in three 
main types (A, B, C), plus a couple of sub-types, 
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Fig. 15. Light curve and color-color diagram for obser- 
vation 40124-01-38-01. Only three PCUs were on (MJD: 
51485.075). 

obtaining a coherent scenario which can be com- 
pared to that of other systems such as XTE J1550- 
564, for which the "ABC" classification was in- 
troduced flWijnands et al. 1999| IRemillard et al. 2002bjl . 
and GX 339-4, where a type-B QPO was found 
flNespoli et al. 2003| ). We now discuss our results in terms 
of relations between derived quantities. 

— Energy dependence: In Fig. 1161 we show the fractional 
rms of the QPO as a function of energy for type-A, 
B and C spectra. The rms amplitude of the QPOs in- 
creases with energy and then flattens above ~10 keV. 
The similarity of the trend in the three QPO types, 
particularly in the B and C cases, is apparent. The 
energy spectra of black hole candidates are often de- 
scribed in terms of two components, one soft/low en- 
ergy component and another hard/high energy com- 
ponent (see e.sr. lTanaka fc Lewin 1995J1 . The observed 
energy dependence of the fractional rms clearly leads 
to the conclusion that the QPOs are associated with 
the high-energy component. Notice that a scenario in 
which the intrinsic fractional rms of the QPOs (i.e. 
the rms normalised only to the source counts from 
the high energy component) is constant - with the ob- 
served decrease in the rms towards lower energy be- 
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Fig. 16. Fractional rms amplitude of the fundamental 
QPO peak as a function of energy for QPO types A (stars, 
obs: 40124-01-12-00), B (40122-01-01-03) and C (40124- 
01-21-00). 



ing caused by a non-variable soft component - can- 
not be ruled out. If this hypothesis were true, the 
QPOs could then be interpreted as bare flux oscilla- 
tions of the high energy component only, without pho- 
ton index changes. Nevertheless, in GRS 1915+105 a 
similar energy dependence for the type-C QPO was 
observed, but the analysis of higher-energy HEXTE 
data showed that the QPO rms decreases above 20 
keV (|Tomsick &: Kaaret 200 fl . Rodriguez et al. (2004) 
found however that this cut-off was not always present, 
and proposed an cxplaination in terms of a contribu- 
tion from the jets to the hard X-ray component. If 
this contribution is significant, and the QPOs are un- 
related to the jet, its presence would thus produce the 
observed cut-off in the energy dependence of the QPO 
rms. 

Time lags: In Figure IT71 we show the behaviour of the 
time lags of the fundamental QPO peak as a function 
of the total rms of the PDS. A separation among the 
three different QPO types is evident: type-C QPOs 
show a strong and well defined trend, with smaller soft 
lags for larger rms, while type-B QPOs have hard lags 
and are clustered around rms of a few percent. The 
latter show also evidence for a dependence of the time- 
lags on the rms. Finally, type-A lags are negative and 
show large scatter. The separation of the three types 
is evident also in Figure IT%1 (left panel), where we plot 
the time lags of the fundamental and of the second 
harmonic. Type-B and type-C QPOs show a similar 
behaviour: the lags of the fundamental and those of the 
second harmonic, with the exeption of a few cases with 
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Fig. 17. Total fractional rms (0.03-64 Hz) vs. time-lags of 
the fundamental QPO peak. 



large error bars, always have opposite signs, resulting 
in a clustering of the points in the second and fourth 
quadrant. This analysis could not be extended to type- 
A QPOs, since they don't show any harmonic peak. 
This peculiar pattern in the time lags of the dif- 
ferent harmonic components was already noticed by 
Remillard et al. (2002b) in the PDS of XTE J1550- 
564. Different lags in these components are difficult to 
interpret. It is of course possible that the oscillations 
at different harmonics have a different physical ori- 
gin, while connected to the same underlying 'clock', in 
which case the difference in sign of time lags would be 
naturally associated to these different physical mecha- 
nisms. However, if one assumes the more natural sce- 
nario in which the harmonic content of the QPO is 
simply the result of the non-sinuoidal shape of the os- 
cillation, the higher harmonics do no have a physical 
meaning by themselves (for a discussion of the sub- 
harmonic see below). Such strong differences in time 
lags would then indicate that the shape of the oscilla- 
tion is different at different energies, in a way which is 
highly reproducible. In order to understand this phe- 
nomenon in detail, a theoretical model for the shape 
of the oscillation is needed. 
— Subharmonic peaks: The right panel of Figure IT%1 giv- 
ing the time lags of the subharmonic versus those of the 
fundamental, shows that subharmonic lags are always 
negative. This behaviour is independent of the QPO 
type and thus of the sign of time lags at the fundamen- 
tal and second harmonic peak frequencies. This some- 
how makes the subharmonic oscillations stand out, and 
suggests a common origin for it in both B and C QPO 
types. 
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Fig. 18. Time-lags of the fundamental QPO peak vs. time-lags of its second harmonic (left panel) and subharmonic 
(right panel). 



The existence of a subharmonic peak opens then 
a serious issue about the physical mechanism that 
would produce them. Even though there are other 
astrophysical examples of subharmonics (see e.g. 
lAikawa fc Antonello 2 000 for the case of Cepheids and 
|Masser &: Tagger 1997| for the case of disk galaxies), 
there are still no generally-accepted physical explana- 
tions. A basic mathematical explaination involves a 
subharmonic resonance, which can arise in forced non- 
linear oscillators. It can be shown (see for example 
Butikov 2002) that for such an oscillator a subhar- 
monic resonance of order n can occur when the driv- 
ing frequency is close to an integer multiple n of the 
natural fundamental frequency. In this scenario, the 
frequency of ~6 Hz seems to play an important role: 
from Figure ^2 it is clear that the coherence of the 
type-B subharmonic peak is higher when the funda- 
mental peak is close to 6 Hz. It is worth noticing how- 
ever that type-C QPOs do not show any peculiar be- 
havior when at this frequency. In a forthcoming paper 
we shall discuss in greater detail the role of this fre- 
quency in different classes of sources, BHCs (see e.g. 
INespoli et al. 2003), Z-sources and atoll-sources (see 
e.g. |Belloni, Parolin fc Casella 2004| >. 



There is of course an alternative scenario in which the 
frequency of the subharmonic peak is in reality the 
fundamental frequency of the system. This hypothesis 
would find some support in the fact that for all types 
of QPOs the subharmonic always shows the same sign 
of lags. If this is the case, there would be no need 
for an explanation of subharmonic peaks; however, 
the subharmonic is not always observed and almost 
never (with the exception of the 'cathedral' cases) the 
strongst peak. Moreover, any model under this sce- 



nario would have to explain the absence of the third 
harmonic. 

Color evolution: In Figure 1191 we plot the Hardness- 
Intensity diagram of the outburst, which shows 
clear analogies with diagrams of other BHCs (see 
IHoman et al. 20011 IBelloni 20031 IR.ossi et al. 2001. 
The energy spectrum is hard at the beginning of 
the outburst, quickly softening as the count rate ap- 
proaches the outburst peak. Here, the flux shows 
strong variations, hardening and softening on short 
time scales. It is during this phase, at the highest count 
rates, that the QPOs were nearly always present. In 
particular, from the right panel of Figure EH it can be 
seen that type-A QPOs appear when the hardness in- 
dicates that the energy spectrum is softer than when 
type-C QPOs are present, while type-B QPOs appear 
at hardness values overlapping those of type-C and 
type-A QPOs, but always restricted between 0.2-0.3. A 
similar correlation between QPO types and hardness 
was found in XTE J1550-564 IjHoman et al. 2001|l . 
From the result reported in Section 4, it is thus clear 
that the type-B QPO can be associated with spec- 
tral transitions. This result is also found in GX 339-4 
(Belloni et al. 2004, in prep.). Either this oscillation is 
excited in a small range of disk parameters (accretion 
rate?), or its frequency depends only weakly on these 
parameters. In either cases, this oscillation seems to 
be an important ingredient in the evolution of the out- 
bursts of a number of BHT. 

After this phase, the count rate first decreased while 
the spectrum still softens, then it showed a hard flar- 
ing episode, after which it continued to decrease while 
maintaining roughly a constant hardness. Finally, af- 
ter more than four months from the beginning of the 
outburst, the energy spectrum hardened, moving to- 
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Fig. 19. Left panel: Hardness-Intensity diagram (see Figure 1 for variables definition). The time resolution is 16 s except 
for the first and the last six observations, which have one point for each observation in order to improve the statistics. 
Different grayscales indicate different QPO types. (A color image is available on line) Right panel: enlargement of the 
region of the diagram where most of type A-B-C QPOs appear. 



wards hardness values comparable to those of the first 
observation. Following Homan & Bclloni (in prep.), 
the roughly horizontal branches at the top and at the 
bottom of Fig. ^j] would therefore correspond to the 
very high/intermediate state (see lHoman et al~ 001). 
appearing at different flux values, while the vertical 
branch to the left would represent the high/soft State 
(as can be inferred also from the low rms values) . Since 
the early part of the rise of the outburst has not been 
observed we can only argue, in analogy with other 
sources, that the source moved along a roughly ver- 
tical branch on the right, which would correspond to 
the hard state. 

6. Conclusions 

From the large set of observations presented here, the 
classification of low-frequency QPOs into three types 
emerges strenghtened and is extended to another source, 
XTE J1859+226. In particular, the type B and C QPOs 
seem to be a key ingredient that is found in a num- 
ber of sources. Their different properties (while at the 
same centroid frequency of ^6 Hz) might provide clues 
to the understanding of the physical mechanisms during 
the VHS/IS of BHCs and are probably related to the two 
"flavors" of VHS/IS observed in these systems (see e.g. 
|Miyamoto et al. ~1 993). The observed energy dependence, 
that rules out a direct disk origin, is an important con- 
strain for the future physical identification of these fea- 
tures. The difference in time-lag behaviour of the different 
peaks and the presence of the subharmonic peak are chal- 
lenging features of the QPO phenomenon. In particular 



our results provide additional evidence that the frequency 
of ^6 Hz is related to some fundamental process. The 
evidence of a treshold triggering the ~6 Hz QPO, super- 
imposed to the roughly exponential decay of the outburst, 
furthermore suggests the existence of a second parameter 
in addition to the mass accretion rate probably responsi- 
ble for the long time-scale evolution of the source. 

In addition, after the analysis of yet another transient 
source, it is clear that the outbursts of these systems, 
in their general evolution, show strong similarities. This 
should be examined in the light of theoretical models for 
the high-energy emission and possibly even used as an- 
other identification tool for these sources. 
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